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PREFACE 


The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under 
Program  Element  62302F,  at  the  request  of  the  Air  Force  Rocket  Propul¬ 
sion  Laboratory  (AFRPL),  AFSC.  The  AFRPL  project  monitor  was  Lt. 

I.  L.  Witbracht,  The  results  were  obtained  by  ARO,  Inc.  (a  subsidiary 
of  Sverdrup  &  Parcel  and  Associates,  Inc. ),  contract  operator  of  AEDC, 
AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project  Number 
V32S-51A.  The  authors  of  this  report  were  W.  D.  Williams  and  J.  W.  L. 
Lewis,  ARO,  Inc.  The  manuscript  <ARO  Ccwtrdl  No.  ARO-VKF-TR-75- 
104)  was  Submitted  for  publication  on  June  30,  1975. 
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1.0  INTRODUCTION  AND  THEORY 


1.1  BACKGROUND 

The  number  of  previous  studies  6f  CO2  condensation  phenomena  in 
gas  dynamic  expansion  flow  fields  is  indicative  of  the  variety  of  applica¬ 
tions  requiring  an  understanding  of  the  process.  Among  such  applications, 
one  finds  the  expansion  of  combustion  products  which  frequently  include 
CO2  and,  also,  the  production  of  high-speed  flow  fields  of  air  which 
naturally  contains  CO2  as  a  significant  constituent.  Previous  studies 
have  utilized  a  wide  variety  erf  flow  diagnostic  techniques  which  have  been 
selected  in  some  cases  because  of  the  intended  application  of  the  results. 

As  an  example,  Hagena  and  Obert  (Ref.  1)  have  employed  a  mass 
spectrometric-retarding  potential  technique  for  studying  the  cluster  size 
distribution  of  the  far -field  region  of  CO2  expansions.  Further  mass 
spectrometric  studies  are  exemplified  by  the  work  of  Bailey  (Ref.  2),  in 
which  not  only  are  the  species  mass  selected  but  also  flow  speed  meas-  ■ 
urements  are  performed  enabling  one  to  obtain  estimates  of  condensate 
mass  fraction.  However,  once  again  it  is  the  expansion  far-field  region 
which  is  studied,  for  at  higher  flow  field  densities  it  is  recognized  that 
skimmer  interaction  effects  do  exist  (Ref.  2)  and  that  the  existence  of 
a  mild  shock  at  the  skimmer  may  be  sufficient  to  dissociate  the  rather 
fragile  clusters. 

Beylich  (Ref.  3)  has  used  the  electron  beam  fluorescence  technique 
to  study  the  density  and  rotational  temperature  of  CO2  plume  expansions, 
but  the  application  was  particularly  tortuous,  for  the  CO2  fluorescence 
spectra  were  not  fully  resolved  and  rotational  temperatures  are  deter¬ 
mined  using  the  emission  band  profile.  Further,  in  the  region  of  conden¬ 
sation  onset  and  growth,  the  gas  density  is  sufficiently  large  to  require 
significant  collision  quenching  corrections. 

Additional  studies  of  CO2  cluster  properties  have  been  performed 
most  notably  by  French  workers  using  the  electron  diffraction  technique, 
and  the  work  of  Audit  (Ref.  4)  is  typical  of  these  efforts.  The  fundamental 
data  of  clusters  provided  hy  this  technique  include  equilibrium,  inter - 
molecular  bonding  distances  of  the  cluster,  and  possibly  the  cluster 
temperature.  However,  the  diagnostic  technique  taken  alone  is  lacking 
in  the  basic  gas  dynamic  data  which  are  necessary  to  unambiguously 
determine  the  cluster  size.  Further,  the  technique  is  a  far-field, 
low -density  experimental  tool  which  suffers  from  the  same  disadvan¬ 
tages  as  does  mass  spectroscopy. 
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Finally^  Rayleigh  scattering  has  been  successfully  empLc^ed  by 
Beylich  (Ref.  5)  and  Lewis  and  Williams  (Refs.  6  and  7)  to  observe  the 
condensation  onset  and  growth  regions  of  a  gas  dynamic  expansion  process. 
The  work  of  Refs.  6  and  7  extended  the  use  of  the  technique  to  include 
observations  of  scattering  depolarization  factors  and  thereby  deduce 
informative  data  regarding  cluster  asymmetry  characteristics.  It  should 
also  be  noted  that  the  work  of  Refs.  6  and  7  approached  the  scattering 
process  from  a  molecular  rather  than  macroscopic  basis  as  did  Beylich, 
and  this  approach  has  enabled  fundamental  information  to  be  obtained 
regarding  the  accretion  process  of  cluster  formation. 

More  recently,  the  work  of  Lewis  and  WiUi^uns  (Ref.  8)  has  demon¬ 
strated  the  application  of  Raman  scattering  to  condensing  flow  field  studies, 
thereby  providing  a  means  of  measviring  local  values  of  both  gas  density 
and  temperature.  It  is  to  be  noted  that,  contrary  to  the  techniques  of 
Refs.  1  through  4,  both  Rayleigh  and  Raman  scattering  are  ap^Aicable 
to  the  high -density  regions  of  the  flow  field  which  are  most  convenient 
for  condensation  studies. 

Despite  the  previous  interest  in  CO2  as  described  in  Refs.  1  through 
5  and  the  references  quoted  therein,  little  direct  information  exists  re¬ 
garding  the  spatial  location  within  the  flow  field  of  condensation  onset,  or 
the  degree  of  supersaturation  capable  of  CO2  and,  also,  the  rate  of  growth 
of  the  clusters  following  condensation  initiation.  Further,  general  scaling 
laws  for  various  types  of  gas  sources  such  as  sonic  orifices  and  conical 
nozzles  have  not  been  verified.  It  should  be  noted  that  CO2  presents  a 
particularly  challenging  case:  not  only  must  condensation  be  considered, 
but  also,  for  very  modest  reservoir  temperatures,  sufficient  vibrational 
mode  excitation  exists  to  require  consideration  of  vibrational  relaxation 
phenomena  as  well.  Consequently,  prediction  of  the  uncondensed  gas 
dynamic  expansion  of  CO2,  which  determines  the  detailed  behavior  of  the 
condensation  process  at  onset,  is  not  a  trivial  problem. 

The  data  of  this  report  present  the  results  of  a  study  of  Rayleigh 
scattering  of  the  CO2  monomers  and  clusters  produced  by  a  sonic  orifice 
expansion  from  ambient  temperature  reservoir  conditions.  Because  of 
the  complexity  of  the  study,  further  work  is  both  in  progress  and  planned, 
and  it  is  hoped  that  the  results  and  analysis  of  CO2  expansions  will  be 
presented  in  more  completeness  in  following  reports. 
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1.2  RAYLEIGH  SCATTERING 

The  basic  equations  and  rationale  for  the  application  of  Rayleigh 
scattering  to  the  study  of  condensing  gas  flow  fields  have  been  presented 
and  discussed  in  previous  publications  (Refs.  .6  and  7).  Consequently, 
only  a  summary  of  the  pertinent  information  and  equations  will  be  given 
in  this  section. 

It  will  be  recalled  from  Refs.  6  and  7  that  for  an  incident  polarized 
laser  source  of  wavelength  (A)  and  intensity  (1^)  interacting  with  a  gas 
sample  of  number  density  (N^)  the  Rayleigh  scattered  intensity  (Ig)  can 
be  written  as 


where  is  the  molecular  polarizability,  N  is  the  gas  reservoir  density 
and  is  a  constant  including  optical,  geometrical,  and  fundamental 
atomic  constants.  If  the  flow  field  is  isentropic,  obviously  Ig/l^  is  a 
direct  measiire  of  the  isentropic  density  ratio  (Nj/Ng)^,  where  the  super¬ 
script  zero  denotes  isentropic  conditions.  As  shown  in  Refs.  6  and  7 
for  a  condensing  flow  field  consisting  not  only  of  monomers  (denoted 
the  subscript  one)  but  also  clusters,  or  i-mers  (denoted  by  subscript  i), 
the  scattered  intensity  can  be  written  as 

/!„)“]  -  (Ki/Nj)  =  (N^/Np  •  .  £  p(i)  .  (a/ai)®  (2) 

i=2 

where  N'p  is  the  local  total  number  density,  is  the  condensate  mole 
fraction  and  p(i)  is  the  probability  distribution  function  of  the  i-mer.  If 
one  approximates  N^/N?  andN-f/Nj  by  unity  and  then  defines  the  left- 
hand  side  of  Eq.  (2)  as  ^e  scattering  function  f,  one  finds 

f  ^  Xg  p(i)(a[/aj)^  (3) 

Further,  assuming  weak  van  der  Waal's  binding  of  the  clusters,  one 
can  approximate  as 

“i  =«  '“1  (4) 

i.  e. ,  cluster  polarizability  is  assumed  to  be  additive.  Then 

f  »  p(i)  (5) 


7 


AEDC-Tn-75-146 


For  a  monodisperse  distribution, 

p(i)  =  1  i  =  J  . 

=  0  i  #  J 

Then. 

<  “  X.  •  (6) 

SO,  if  knowledge  of  the  mole  fraction  can  be  obtained,  cluster  size  J 
follows. 

Obviously,  prior  to  condensation  onset,  f  is  zero  identically,  and 
positive  deviations  of  f  from  zero  are  Indicative  of  the  initiation  of 
clustering.  Further,  the  parameter  J  of  Eq.  (6)  will  be  a  space -dependent 
quantity  throughout  the  condensation  growth  region,  as  will  X^,  and  f 
is  an  ambiguous  measure  of  the  increase  in  the  mean  cluster  size  and 
mole  fraction  of  condensate. 

The  results  presented  in  this  report  will  be  the  variation  of  the 
scattering  function  f  with  various  experimental  parameters.  The  conden¬ 
sation  onset  and  growth  regions  will  be  clearly  shown  as  will  an  impor¬ 
tant  orifice  diameter/ reservoir  pressure  scaling  law. 


2.0  EXPERIMENTAL  APPARATUS  AND  METHOD 
2.1  GAS  SOURCE  AND  VACUUM  CHAMBER 

The  sonic  orifice  sources  were  of  1.  325-  and  3.  2 -mm  diameter 
with  a  diameter-to-thidcness  ratio  greater  than  20.  The  orifice  was 
located  at  the  end  of  a  tubular,  stainless  -steel  reservoir  of  approxi¬ 
mately  1-cm  inner  diameter,  and  equilibration  of  the  gas  and  tube 
wall  temperatures  was  ensured.  Reservoir  temperature  and  pressure 
were  measured  using  calibrated  transducers,  and  the  flow -field  data 
were  acquired  using  Coleman  grade  CO2.  Two  0.025-ium  particulate 
filters  were  located  in  the  inlet  gas  line  to  minimize  heterogeneous 
condensation  processes.  The  gas  source  was  mounted  on  a  traversing 
mechanism  to  provide  three  degrees  of  freedom  and  to  enable  flow 
field  measurements  to  be  made  using  fixed  optical  Instrumentation. 
The  movement  mechanism  has  an  accuracy  and  reproducibility  of 
0.  013  cm  in  the  axial  direction. 
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The  4-  by  10-ft  Research  Vacuum  Chamber  enclosing  the  motor- 
driven  traversing  mechanism,  as  shown  in  Fig.  1,  was  cryopumped 
using  liquid  N2  and  20  K  gaseous  He.  Background  pressures  of  less 
than  10“^  torr  were  maintained  for  the  maximum  flow  rate  investigated. 
Pressure  and  temperature  fluctuations  in  .the  reservoir  were  negligible. 
A  more  complete  description  is  to  be  found  in  Ref.  6. 
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Figure  1.  Schematic  of  4-  by  10-ft  Research  Chamber. 


2.2  RAYLEIGH  SCATTERING  APPARATUS 

The  optical  and  electronic  configurations  of  the  experimental  appa¬ 
ratus  are  shown  in  Fig.  2,  and  the  polarization  vector  orientation  for 
this  scattering  study  is  shown  in  Fig.  3.  An  argon  ion  laser  of  1.0  w 
intensity  at  514.  5  nm  was  empiloyed  as  the  incident  radiation  source. 

The  collection  optics  train  included,  in  addition  to  collection  and  focusing 
lenses,  a  half-wave  plate  polarization  rotator,  and  HN-22  Polaroid®for 
definition  of  the  state  of  polarization  of  the  scattered  radiation.  The 
system,  which  is  described  in  detail  in  Ref.  8,  utilized  a  Spex  1402 
double  spectrometer  with  two  ganged  102-  by  102 -mm  gratings  of  1200 
grooves /mm  and  blazed  at  500.  0  nm.  The  focal  length  of  the  instrument 
was  0.  85  m  with  a  reciprocal  linear  dispersion  of  approximately  0.  44 
nm/mm.  The  three  slit  widths  were  set  at  200,  400,  and  400  pm  in 
order  of  input  to  output  slits.  Photon  counting  electronics  in  conjunction 
with  a  thermoelectrically  cooled  photomultiplier  tube  served  for  data 
acquisition.  The  experimental  procedures  and  method  used  are  presented 
in  Refs.  7  and  8. 
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Figure  3.  Polarization  direction  diagram. 
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3.0  RESULTS  AND  DISCUSSION 


Sonic  orifice  expansions  of  CO2  were  produced  using  source  diam¬ 
eters  (D)  of  1.325  and  3.2  mm.  The  reservoir  pressure  (Pq)  was  varied 
over  the  range  of  50  to  350  torr  for  the  3.  2-mm-diam  source,  and  for 
the  1.  325-mm-diam  source,  the  Pq  range  was  77.  7  to  543,  9  torr.  The 
values  of  P^  for  the  smaller  diameter  source  were  selected  to  achieve 
constant  values  of  p2d  for  the  two  sources  at  e^ch  value  of  Pq  studied. 
The  axial  variation  oi  the  Rayleigh  scattered  intensity  was  determined 
over  the  axial  distance  range  of  0.  8  ^  x/D  ^  60  where  x  is  the  axial 
distance  measured  downstream  from  the  flat-faced  sonic  orifice  source. 
All  cases  studied  were  for  ambient  CO2  reservoir  temperature,  which 
was  nominally  280  K.  Shown  in  Figs.  4  through  8  are  the  experimental 
Rayleigh  scattering  results  for  this  study,  and  the  ordinate  is  (Ref.  6): 


.  (N. 


,/N.) 


i=2  *  ^  >  o 


(7) 


where  is  determined  experimentally  by  an  in  situ  calibration. 

It  is  seen  from  Eq,  (7)  that  in  the  event  of  no  condensation  Ig/l^  will 
vary  directly  as  Nj/N^  which,  in  the  absence  of  vibrational  relaxation 
effects,  is  the  isentropic  variation  (N^/Nq)®.  Also  shown  in  Figs.  4 
through  8  is  the  isentropic  variation  of  Nj/N^  as  given  by  the  Ashkenas- 
Sherman  (Ref.  9)  equations  for  7  =  1.4,  The  three-term  expansion  of 
Ref.  9  was  used  to  enable  accurate  calculations  for  the  axial  distance 
region  of  x/D  >  1.  0.  The  axial  locations  at  which  saturation  occurs 
[  (x/D)a  ]  were  estimated  for  CO2  using  the  data  of  Ref.  10  and  as¬ 
suming  an  isentropic  expansion  of  CO2  characterized  by  a  specific 
heat  ratio  (7)  of  1.  40.  Differences  of  approximately  ±20  percent  in 
(x/D)g  result  from  a  variation  in  7  of  approximately  ±5  percent.  More 
accurate  values  of  (x/D)g  are  not  required  for  this  report,  and  it  was 
found  that  for  all  cases  studied 


(X/D),  0.4  (8) 

Consequently,  the  rate  equations  which  describe  the  onset  and  growth 
regions  of  CO2  must  be  applicable  in  the  near  field  of  the  expansion. 

It  is  seen  from  Fig.  4  that  the  expansion  with  the  lowest  Pq  yielded 
agreement  with  the  7  =  1. 40  isentropic  prediction  of  Ref.  9  for  x/D  ^3. 
Further,  from  Figs.  4  through  8,  one  sees  that  an  increase  in  Pq  resulted 
in  condensation  onset  and  that  the  onset  location  decreased  while  the  mag¬ 
nitude  of  deviation  of  1q/Iq  from  the  isentropic  dependence  increased  as 
pQ  increased. 
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Figure  4.  Axial  variation  of  scattered  intensity  (U/lo)  at 
Pq  =  50  and  77,7  torr. 
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Figure  6.  Axial  variation  of  scattered  intensity  (It/lo)  at 
Po  =  200  and  310.8  terr. 
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Additionally,  one  notes  that,  although  the  data  from  the  two  sonic 
orifice  sources  are  in  approximate  agreement  for  constant  values  of 

there  is  an  observatile  discrepancy  in  these  data.  Consequently,  P^D 
scaling  yields  an  approximate  scaling  law  for  CO2,  but  further  refinements 
are  required  to  the  scaling  law  to  reduce  the  existing  discrepancy. 

The  scattering  function  flp^)  where  the  superscript  and  subscript 
denote  the  orifice  diameter  and  reservoir  pressure,  respectively,  is 
shown  in  Fig.  9  as  a  function  of  for  the  3.  2'mm-diam  source.  The 
rapid  increase  in  f  following  onset  is  quite  obvious  as  are  the  orders  of 
magnitude  increase  in  f  as  the  reservoir  pressure  increases.  The  quan¬ 
titative  definition  of  condensation  onset  can  take  several  forms.  Used  in 
this  work  are  two  criteria  which  involve  either  the  scattering  function  f 
or  the  scattered  intensity  Clg/Ip).  The  first  method  employs  the  criterion 
of  a  ten  percent  increase  of  Is/Iq  ^^l^'tive  to  its  isentropic  value,  and  all 
condensation  parameters  derived  using  this  criterion  are  denoted  by  the 
subscript  0  .  The  second  method  uses  the  f -function;  and  specifically, 
the  variation  of  f  with  x/D  is  used  to  extrapolate  the  f  function  to  the 
onset  location.  The  value  of  the  condensation  parameters  appropriate 
to  this  extrapolated  axial  position  at  which  f  equals  zero  are  denoted  hy 
the  subscript  9 .  Then  the  super  saturation  parameters  S  and  S'  are 
defined  as 

-  T,  -  (9) 

and 

%4.  - 


where  the  subscript  oodenotes  free-stream  values.  Figures  10  and  11 
show  the  dependence  of  the  supersaturation  parameters  on  reservoir 
pressures  for  the  two  sonic  orifice  sources.  From  Figs.  10  and  11,  it 
is  seen  that  the  temperature  supercooling  values  determined  from  Sq 
and  5^  are  not  substantially  different  and  are  on  the  order  of  100  K.  The 
locus  of  the  saturation  temperature  Tg  with  is  shown  in  Fig.  10. 
y  The  variation  of  the  pressure  supersaturation  ratios  is  seen  to  decrease 
significantly  as  increases. 

Finally,  empirical  functional  relations  of  f  with  axial  distance  X/D, 
(X/D}g  ,  Pq,  and  D  were  obtained  using  the  results  of  the  axial  variation 
of  the  scattering  function  data. 
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Figure  9.  Axial  variation  of  scattering  function  f  jpjj 
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Figure  11.  Variation  of  supersaturation  paramatars  (S^  and  S^)  with  Pq. 

The  result  found  for  the  cases  studied  in  the  onset  region  is 

=  ben  [{x/d)/(x/d)(9] 

where 

b  =  (12) 

and  Cq  is  a  constant.  After  cessation  of  growth,  the  scattering  function 
was  found  to  be  described  by 

^  (13) 

where  is  a  constant. 

4.0  CONCLUSIONS 

The  results  of  Rayleigh  scattering  measurements  of  CO2  sonic 
orifice  expansions  have  demonstrated  the  onset  of  condensation  and 
subsequent  cluster  growth  as  well  as  the  very  large  increase  in  the 
scattering  function  (f)  with  increasing  reservoir  pressure.  The  results 
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for  very  low  Pq  and  un  condensed  CO2  indicate  that  the  far -field  expan¬ 
sion  region  is  described  by  the  isentropic,  7  =  1.4  predictions  of  Ref.  9, 
However,  for  all  cases  studied,  the  saturation,  onset,  and  growth  regions 
for  condensation  all  occur  in  the  more  complicated  region  of  the  expan¬ 
sion  where  the  flow -field  is  yet  developing,  not  in  the  far-field  region 
of  the  expansion  where  one  may  use  various  asymptotic  relations  for  the 
gas  dynamic  parameters.  Empirical  relations  of  the  scattering  function 
dependence  on  Pq,  D,  and  X/D  were  found  and  will  serve  as  a  basis  for 
predicting  gross  condensation  characteristics. 


REFERENCES 

1.  Hagena,  O.  F.  and  Obert,  W,  "Cluster  Formation  in  Expanding 

Supersonic  Jets:  Effect  of  Pressure,  Temperature,  Nozzle 
Size,  and  Test  Gas,  "  The  Journal  of  Chemical  Physics, 

Vol.  56,  No.  5,  March  1972,  pp,  1793-1802. 

2.  Bailey,  A.  B.  "Effects  of  Condensation  of  Gas  Velocity  in 

a  Fr ee- Jet  Expansion.  "  AEDC-TR-73-93  (AD762503), 

June  1973. 

3.  Beylich,  A.  E.  "Experimental  Investigation  of  Carbon  Dioxide 

Jet  Plumes. "  The  Physics  of  Fluid,  Vol.  14,  May  1971, 
pp.  898-905. 

4.  Audit,  P.  "Liasons  Intermolecularies  dan  les  Jets  Supersoniques 

Etude  par  Diffraction  d'Electrons,  "  Journal  Physique,  Vol.  30, 
1969,  pp.  192-202. 

5.  Beylich,  A.  E.  "Condensation  in  Carbon  Dioxide  Jet  Plumes.  " 

AIAA  Journal,  Vol.  8,  May  1970,  pp,  965-967. 

6.  Lewis,  J.  W.  L. ,  and  Williams,  W.  D.  "Argon  Condensation  in 

Free -Jet  Expansions."  AEDC-TR-74-32  {AD782445),  July  1974. 

7.  Lewis,  J.  W.  L. ,  Williams,  W.  D. ,  Price,  L.  L. ,  and  Powell, 

H.  M.  "Nitrogen  Condensation  in  a  Sonic  Orifice  Expansion 
Flow."  AEDC-TR-74-36  (AD783254),  July  1974, 

8.  Lewis,  J,  W.  L. ,  and  Williams,  W.  D,  "Measurement  of 

Temperature  and  Number  Density  in  Hypersonic  Flow  Fields 
Using  Laser  Raman  Spectroscopy.  "  AIAA  Paper  75-175, 

AIAA  13th  Aerospace  Sciences  Meeting,  Pasadena,  California, 
Januaiy  1975. 


21 


AEDC-TR-7B-146 


9,  Ashkenas,  H. ,  and  Sherman,  F.  S.  "The  Structure  and  Utilization 
of  Supersonic  Free  Jets  in  Low  Density  Wind  Tunnels.  " 
Rarefied  Gas  E^namics,  Edited  by  J.  H.  de  Leeuw,  Fourth 
.  Symposium,  Vol.  n, ,  Academic  Press,  New  York,  1966, 
pp.  84-105. 

10.  HUsenrath,  J. ,  Beckett,  C.  W. ,  Benedict,  W.  S, ,  Fano,  L. , 

Hoge,  H.  J.,  Masi,  J.  F. ,  Nuttall,  R.  L. ,  Touloukian,  Y.  S. , 
and  Woolley,  H.  W.  Tables  of  Thermal  Properties  of  Gases. 
National  Bureau  of  Standards  Circular  564,  Washington^  D.  C. , 
November  1955, 

NOMENCLATURE 


D  Orifice  diameter 

f  Scattering  function 

li^  Incident  laser  intensity 

Ig  Scattered  intensity 

Ka  Optical  constant 

N  Number  density 

N^  Number  density  of  i-mer 

P  Presstire 

p(i)  Probability  distribution  function  of  i-mer 

S  Supers aturation  parameters  defined  Eqs.  (9)  and  (10) 

T  Temperature 

Mole  fraction  of  condensate 
X,  y,  z  Directional  coordinates 

Polarizability  of  i-mer 
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y  Specific  heat  ratio 

X  Wavelength 

SUBSCRIPTS 

c  Condensate 

s  Saturation  value 

T  Total  value 

o  Reservoir  parameter 

0  Defined  by  Eqs.  <9)  and  {10} 

Defined  by  Eqs.  (9)  and  (10) 

no  Free -stream  value 

SUPERSCRIPT 

0  Isentropic  value 

Note;  Unless  otherwise  noted,  cgs  units  are  used. 
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